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Abstract. We present an analysis of an XMM-Newton observation of the M dwarf binary EQ Pegasi with a 
special focus on the the spatial structure of the X-ray emission and the analysis of light curves. Making use of 
data obtained with EPIC (European Photon Imaging Camera) we were for the first time able to spatially resolve 
the two components in X-rays and to study the light curves of the individual components of the EQ Peg system. 
During the observation a series of moderate flares was detected, where it was possible to identify the respective 
flaring component. 



1. Introduction 

X-ray observations with the Einstein Observatory and 
ROSAT have shown the ubiquitous occurrence of coronae 
around most classes of stars. ROSAT studies of volume- 
limited complete samples of cool stars in the immediate 
solar neighborhood have shown coronal formation around 
late-type cool dwarf stars with outer convection zones to 
be universal; all stars investigated with sufficient sensi- 
tivity we re found to b e surrounded by X-ray emittin g 
coronae l)Schmitt et all Il995l ISchmitt fc LiefkeL 12003^ . 
Interestingly, fully convective M dwarfs have also been 
found to be very active with frequent flares. 

EQ Peg is a nearby (6.25 pc) visual binary (period 
~ 180 yr, separation 5.2") consisting of two M dwarfs of 
spectral type M3.5 and M4.5. It was f irst o bserved pho- 
toelec tricallv to flare bv iRoauesl l)l954|) . and lOwen et alJ 
found both components of the system to be flare 

stars. 

EQ Peg has been observed at radio, optical, EUV, and 
X-ray wavelengths. Observations in the optical focused on 
the flare nature of EQ Peg and marked emission li ne vari- 
abilit y during photometric quiescence was found l|Bot)ri , 
as well as freq uent optical flares on both compo- 
nents l)Rodoncl Il978fl. A VLA m a p of E Q Peg at 6 cm 
was presented by iTopka fc Marshl jl982() . They resolved 
both components and interpreted the radio emission as 
"quiescent" since they found it unlikely that both com- 
ponents flared at the same time. The radio emission was 
confined to each component and ITopka fc Marshl £l982) 
concluded that radio production mechanisms do not de- 
pend on binary interaction (which is plausible due to the 
separation of ~ 25 AU). 
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EQ Peg was observed by all major previous imag- 
ing X-ray missions and again found to flare frequently. 
EQ Peg was obser ved with the Einstein Observatory 
l|Vaiana et allll98ll) and is cont ained in the ESS (Einstein 
Slew Survey) l|Elvis et allll992jh EXOSAT detected an in- 
tense long durat ion flare during a coord inated observation 
with the VLA l|Pallavicini et all 11986^1 ; a detailed mod- 
elling of these fl ares and the underly ing physical properties 
is presented bv lPolotto et all (11988). E Q Peg was also de- 
tected in the ROSAT all-sky survey (|Hunsch et allll999j) 
and rapid flaring was simultaneously observed at opti- 
cal and X-ray wavelen gth with MEKASPEC and ROSAT 
l|Katsova et allEoO^ . where the source brightened in X- 
rays by a factor of ~ 15. A coordinated VLA, optical, 
EU VE, and RXTE mo nitoring of EQ Peg was carried out 
bv iGaene et all l[l998). They found a classic stellar flare 
with a rapid impulsive phase (radio burst) followed by 
rapid chromospheric heating and cooling (U-band) and 
more gradual coronal cooling (X-ray and extreme-UV). 
In addition they found atypical flares with either highly 
polarized emission with no counterparts at shorter wave- 
lengths or moderately polarized flares that often have 
shorter- wavelength counterparts . 

EQ Peg was also observed with XMM-Newton. In 
Sect.|21wc describe the observations and the methods used 
for data analysis. Here we focus on the data from the EPIC 
instruments in order to obtain spatial and temporal infor- 
mation on the two components of EQ Peg. In Sect. |3| we 
present the results followed by a summary and discussion 
in Sect.H 

2. Observation and data analysis 

EQ Peg A/B (V=10.32mag/12.4mag) was observed on 
2000 July 9 (MJD=51734) with XMM-Newton. The 
15ksec observation of EQ Peg (see Tab.^l provided use- 
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ful data in all EPIC (European Photon Imaging Camera) 
detectors. The EPIC instrument consists of three CCD 
cameras with two different types of CCD design, resp. 
two MOS (Metal Oxide Semi-conductor CCDs) and one 
PN (pn CCDs), providing imaging and spectroscopy. The 
EPIC cameras offer the possibility to perform extremely 
sensitive imaging observations over the telescope's field of 
view of 30 ' and in the energy range from 0.15 to 15 keV 
with good angular and moderate spectral resolution. A de- 
taile d descrip t ion of the XMM instruments can be found 
in lEhle et all l|2003|) . All EPIC instruments (MOS/PN) 
operated in the full frame mode with the thick filter in- 
serted. Unless otherwise indicated we used for our pur- 
poses the full energy bandpass of the EPIC instruments, 
resp. 0.15 keV to 12.0/15.0keV. 



Table 1. Observation log of EQ Peg 

Instrument (Mode) Duration (s) Obs-timc 

MOS (FF.thick F.) 14600 2000-07-09Tll:39:13-15:42:32 

PN (FF.thick F.) 12410 2000-07-09T12:20:16-15:47:05 



The data were reduced with the standard XMM- 
Newton Science Analysis System (SAS) software, version 
5.4.1. Light curves and images were produced with stan- 
dard SAS tools and standard selection criteria were ap- 
plied for filtering the data. In Fig. ^ we show the image 
obtained with the MOS1 detector. The image looks elon- 
gated and it is reasonable to assume the elongation is due 
to emission from both components of EQ Peg. The image 
elongation can be seen clearly in MOS1, but not in MOS2 
and PN due to the triangular shape of the point spread 
function (PSF) for MOS2 and the bigger pixel size in PN. 
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Fig. 1. Image of the EQ Peg system (MOS1) with linear 
brightness scaling. The image is elongated, suggesting the 
presence of two components. Analysis shows that the X- 
ray brighter component is EQ Peg A. 



For a quantitative analysis of the MOS1 image we de- 
veloped a fitting procedure applicable to the measured 
event distribution in order to confirm the detection of 
the two components of EQ Peg and to determine accu- 
rate source positions and count rates. For this procedure 
we optimize a set of parameters describing the modelled 
event distribution on the sky plane. The modelle d event 
distri bution is constructed on the basis of the PSF l|Kirscbl 
120021) . which is composed of a King component plus back- 
ground 



PSF = 



A 



(l + (^) 2 ) a 



+ BKG 



(1) 



with the core radius r c , the slope a (model parameters), 
the distance to the peak position r, and the amplitude 
A (source parameters). The background flux BKG was 
determined from source free regions in the detector and 
is kept as a fixed model parameter. We place two such 
King components near the positions of the two sources and 
apply an optimization algorithm seeking best-fit values for 
the amplitudes Ai and A2 and the positions r\ = (x\, y\) 
and T2 = {xi, y-i) of each component. We thus determine 
the model counts in the spatial bin 
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n; i = + + BKG (2) 

- (1 + (^) 2 ) Q (l + (^) 2 )« 

and apply Powell's algorithm ijPress et al.lll9 99). a robust 
multi-dimensional minimizing routine, in order to mini- 
mize the likelihood function 



(3) 
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with the measured counts Cjj and model counts The 
model counts are constrained to conserve the total 
number of counts, i.e., n id = £\ ^ c it j. 

The mean PSF-para meters fry, and a ) as determined 
by in-flight calibration l|Ghizzardi l200l[) did not lead to 
good fit results for our EQ Peg data. It turned out that 
especially for the value of the core radius (r c ) binning and 
pile-up effects have to be considered, while the variation 
of the slope a is only moderate. In order to find a bet- 
ter representation for the shape of the PSF we carried 
out a parameter study for our data from EQ Peg and for 
XMM-Newton observations of single point sources, e.g., 
e Eri with a comparable detector configuration and pile- 
up level. From this study we found a slope of a — 1.45, 
which agrees with the calibration value, and a value for 
the core radius of r c = 4.45, i.e., a reduction of ~ 20%, to 
be better suited to model our data. 

With the redetermined PSF parameters we modelled 
the event distribution. In Fig.[21we show a one-dimensional 
representation of our fit results for the EQ Peg observa- 
tion. For this purpose we integrated along the declination 
axis, which almost matches the main axis of the EQ Peg 
system and binned the data into a histogram. As can be 
seen in Fig. O the model describes the data very well. 
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3.2. Identification of flare activity 




arcsec 

Fig. 2. The event distribution of EQ Peg, fitted with the 
PSF model. Shown are the individual components as well 
as the sum of both components compared to the data (his- 
togram). 



3. Results 

3.1. Determination of exact source positions 

A first inspection of Fig. ^ suggests the presence of two 
sources. Since the angular resolution of the MOS1 detector 
is of the same order as the separation of the two compo- 
nents of EQ Peg (~ 5"), it is clear that we are operating 
at the limit of the spatial resolving power of the MOS1 
detector. We applied our PSF algorithm to the EQ Peg 
dataset to determine precise source positions. We use the 
MOS1 data from a 50x50" field centered on the position 
of the EQ Peg system. The calculated source positions 
as listed in Table |21 agree well within th e errors with the 
optic al positions taken from literature jPerrvman et all 
119971) . Here proper motion corrections were not applied 
since they are small because the observation took place in 
July 2000. We therefore conclude that we indeed identified 
the two X-ray sources with the optical counterparts. 

The algorithm was also used to determine the abso- 
lute number of counts per source and the count ratio of 
the sources. EQ Peg A was on average a factor of ~ 3.5 
brighter than EQ Peg B during the total observation. 



Table 2. Position fit results (la errors), ref. data from 
Simbad (FK5/2000) 





EQ Peg 


A 


EQ Peg 


B 




Fit 


Ref. 


Fit 


Rcf. 


RA (23:31:) 


52.17 ±0.01 


52.18 


52.57 ±0.02 


52.53 




(352.9674°) 




(352.9690°) 




DEC (+19:56:) 


14.10 ±0.05 


14.10 


14.05 ±0.10 


13.90 




(19.9373°) 




(19.9372°) 




EQ Peg System 
Separation (") 


Fit 
6.0 ±0.3 




Ref. 
5.2 





Count ratio A/B 3.4 ± 0.2 4.2 ± 0.3 (quics.) 2.9 ±0.2 (flare) 



From the image obtained with MOS1 the two components 
of EQ Peg can be separated. In Fig. |3| we plot the light 
curves for the EQ Peg system as observed with the differ- 
ent EPIC detectors. Inspection of the total light curves in 
Fig. |3 shows that the EQ Peg system stayed more or less 
quiet during the first 8 ksec, afterwards a rise in count rate 
is detected in all three detectors. We therefore divided the 
data set into two parts separated at t=7.86 ksec where we 
consider the first part the quiescent phase and the second 
part the flaring phase. 



RATE (count/s) 




5 10 15 



TIME (ksec) 

Fig. 3. Light curve of the EQ Peg observation as measured 
by the PN (black) and MOS1+2 (grey) detectors with 
100 s binning. 

In Fig. 01 we display the event distributions of these 
two subsets; the histograms were created in the same way 
as in Fig. |21 and are corrected for the different integration 
times. Here EQ Peg A is the X-ray brighter component 
on the right, EQ Peg B corresponds to the weaker com- 
ponent on the left edge of the event-distribution. In the 
quiescent phase EQ Peg A dominates the emission and the 
second component is only marginally visible, while during 
the flaring phase EQ Peg B brightens up and becomes 
more clearly visible. 

The analysis of the X-ray images as carried out for the 
total observation fSect. l3~T)l can be repeated for the differ- 
ent phases of activity. Application of our PSF algorithm 
to the two subsets with variation of only the amplitude 
parameters returned the count rate of EQ Peg A to be 
a factor > 4 higher than for EQ Peg B in the quiescent 
phase, while in the flaring phase the ratio was < 3. The 
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Fig. 4. Comparison of the event distribution during the 
quiescent and flaring phase. 



results of the fitting procedure for the EQ Peg observation 
are summarized in Table [21 

3.3. Reconstruction of individual light curves 

Having found that the major flaring activity is due to EQ 
Peg B, we decided to carry out a systematic light curve 
analysis of both components of the EQ Peg system. As 
a first approach we extracted two different light curves 
from the MOS1 data by placing a circular region with 
2.5" radius around each component, concentrating on the 
core of the PSF. 



the quiescent emission level. Clearly, the extraction re- 
gions used do contain photon contamination from the re- 
spective other component due to the wings of the PSF. 
Nevertheless, the individual light curves also suggest un- 
corrected variability between the A and B components. In 
particular, the flaring at the end of the observations (e.g, 
t > llksec) seems to originate from the A component. 



RATE (count /s) 
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1 5 



time (ksec) 

Fig. 6. Division of the MOS1 light curve into seven time 
intervals (top) and derived flux ratios (bottom). 




TIME (ksec) 

Fig. 5. Light curves of the two components extracted from 
a circular region around each component. 

In Fig. Owe show these light curves binned into 7.5 min 
bins. Again, EQ Peg B is more variable and its count rate 
rises by a factor of ~ 2 at the peak of the flare while 
EQ Peg A shows only marginal brightening compared to 
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Fig. 7. Light curves of EQ Peg B (top) and A (bottom) 
calculated with the PSF fitting algorithm. 

For a more detailed quantitative treatment we utilize 
our PSF algorithm in order to reconstruct the individual 
light curves. We divide the dataset into seven time inter- 
vals covering the various phases of activity as shown in the 
upper panel of Fig. H3 From our PSF algorithm a count 
ratio for each time interval can be determined and in the 
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bottom panel of Fig. we show the development of this 
count ratio. While the main flaring activity is located on 
EQ Peg B (indicated by the decrease of the A/B count ra- 
tio after ~ 5 ksec) there is also some activity on EQ Peg A 
especially during the later phase of the observation. 

From these ratios we calculated light curves for each in- 
dividual component, which are shown in Fig.0 Although 
the light curves consist of rather large time bins the main 
features visible in Fig.Elare also present, i.e., a flare on EQ 
Peg B around 10 ksec and the flaring activity on EQ Peg A 
towards the end of the observation. The rise in count rate 
associated with the flare activity on both components is 
nearly equally strong, i.e., ~ 0.2counts/s, however, the 
relative change in count rate is much higher on EQ Peg B. 



3.4. Development of the spectral hardness 

In order to quantify eventual changes in the physical con- 
ditions accompanied by rising count rates we calculate 
a spectral hardness ratio for the sum of both compo- 
nents in two energy bands, resp. 0.2 — 0.5 keV (soft) and 
1.0 — 10.0 keV (hard). The hardness ratio was calculated 
from PN data which were cleared for pile-up effects. 

In Fig. |H1 we show the spectral hardness ratio 
(hard/soft) for this observation binned every five min- 
utes vs. the measured count rate. The hardness ratio in- 
creases during the times where we also observe increases 
in count rate. We calculated a linear correlation coefficient 
and found a correlation probability of 99.9%. We there- 
fore interpret the simultaneous increase in count rate and 
hardness of the emission as flare heating of the coronal 
plasma. When comparing the calculated average hardness 
ratios during the time intervals defined in Sect. 13.21 as the 
quiescent phase and the flaring phase we find an increase 
in hardness by 12±2 % from the quiescent phase into the 
flaring phase, whereas the increase during individual flares 
is substantially stronger. 



HR (hard/soft) 
_L 




RATE (counts/s) 

Fig. 8. Hardness ratio (l-10keV/0.2-0.5keV) derived 
from PN data vs. count rate. 



The average X-ray luminosity between 0.2 keV and 
10.0 keV was calculated to be 4.2 x 10 28 ergs/s from ap- 
plication of spectral models. We compared the X-ray lu- 
minosities obtained from a number of different reasonable 
models and find consistent results. In comparison with 
previous measurements the energy flux during this obser- 
vation was ~ 20% ab ove the flux obtained with ROSAT 
l|Hiinsch et all Il999(l and ~ 50% below the values ob- 
tained with Einstein l|Vaiana et "atl. ll98lf) in the respec- 
tive energy bands. 

4. Summary and Discussion 

Our analysis of EQ Peg is another example of how 
the high-resolution X-ray telescopes Chandra and XMM- 
Newton allow to resolve sources down to a un- 
precedented spatial reso lutio n (for other examp les see 
IStelzer fc BurwTtl $mk andlAuda.rd et^hl ffijjjjft). 

With the XMM-Newton observation of the EQ Peg 
system we were able to separate the two components for 
the first time in X-rays. Using a PSF model fit procedure 
we can reconstruct the source positions and show that 
both components are flaring X-ray emitters. On average, 
we found the A component brighter by a factor ~ 3.5 for 
the total observation. 

During this observation a series of medium flares was 
detected. We were able to determine count ratios for EQ 
Peg A/B for the different phases of activity. During the 
early (quiescent) phase of the observation the emission is 
strongly dominated by EQ Peg A, which is a factor of 
^4 — 5 brighter than EQ Peg B. Comparison of the quies- 
cent and active phases made it possible to associate most 
of the flaring with EQ Peg B, which nearly doubled it's 
X-ray brightness during the peak of the flare. The count 
ratio during the peak of the flare on EQ Peg B dropped 
to ~ 2.5 — 3. We also found evidence for flaring activity 
on EQ Peg A towards the end of the observation, consis- 
tent with prev ious find i ngs th at both stars exhibit flaring 
behavior fe.g.. lRodonclll978j) . In fact, the relative bright- 
ening during the flares is much stronger for EQ Peg B, but 
the absolute increase in flux is comparable for both stars. 
The energy released by these flares is obviously very simi- 
lar, although the quiescent emission level is quite different. 
The flaring X-ray emission of the EQ Peg system shows 
the typical hardening in the spectral energy distribution 
as expected for stellar flares. 
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